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We study the critical behavior of Boolean variables on scale-free networks with competing inter- 
actions (Ising spin glasses). Our analytical results for the disorder-network-decay-exponent phase 
diagram are verified using large-scale Monte Carlo simulations. When the probability of positive 
(ferromagnetic) and negative (antiferromagnetic) interactions is the same, the system undergoes a 
finite-temperature spin-glass transition if the exponent that describes the decay of the interaction 
degree in the scale-free graph is strictly larger than 3. However, when the exponent is equal to or less 
than 3, a spin-glass phase is stable for all temperatures. The robustness of both the ferromagnetic 
and spin-glass phases suggest that Boolean decision problems on scale-free networks are quite stable 
to local perturbations. 

PACS numbers: 75.50.Lk, 75.40.Mg, 05.50.-fq, 64.60.-i 



I. INTRODUCTION & MOTIVATION 



Networks play an integral role in all fields of science, 
as well as numerous industrial applications. Virtually 
any process or group of interacting entities can be de- 
scribed by a network. In particular, the study of scale- 
free networks — i.e., networks where the degree distribu- 
tion follows a power law — found a renewed interest in 
the last decades after it had been shown that the in- 
ternet follows such a network topology Since then, 
there have been numerous studies showing that a mul- 
titude of networks ranging from computer networks, to 
protein-protein interaction networks, semantic networks 
and, in particular social networks such as citation net- 
works or sexual partner networks are also well-described 
by scale-free networks Q. 

Scale free networks have edge degrees {ki} distributed 
according to a power law A, with the probability pk for 
a node to have k neighbors satisfying 




pk a 



(1) 



A typical network is shown in Fig. [TJ While few nodes 
have many edges connecting them to other nodes, many 
nodes have few edges - the distribution of these following 
a power law. Although there have been several studies of 
Boolean variables on scale-free networks with social in- 
teraction networks in mind, most studies have focused on 
'friend' networks such as, for example, the facebook Q 
network where person A can 'friend' person B. Friend- 
ship can then be defined via a network edge between A 



FIG. 1: (Color online) Typical network simulated. The con- 
nectivity matrix is selected according to Eq. ((T|, i.e., the edge 
degrees are distributed according to a power law. This ex- 
plains why few nodes (larger circles) have many connections 
(darker lines), while many nodes (small dots) have fewer con- 
nections (lighter lines). Note that the minimum connectivity 
in our simulations is 3, i.e., each node has at least 3 neighbors 
and the maximum connectivity scales with the square root of 
the number of nodes. Data for A'' = 512 nodes. 
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and B with a positive weight. However, other networks 
exist where two persons A and B can either be 'friends' 
or 'foes,' i.e., a network with both positive- and negative- 
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weight bonds. This appears, for example, in the slashdot 
network Q or when studying the robustness of opinion 
formation in, e.g., an election process. The latter type of 
network is rarely studied, possibly due to the difficulties 
introduced by the negative- weighted edges in the system. 
However, they find wide applicability to many fields of 
science such as the aforementioned social networks, as 
well as other applications such as interaction networks 
between proteins or genes. 

Why are Boolean problems on these scale- free networks 
interesting? Because they can be seen as the simplest 
model to study how general consensus forms on such a 
network for a decision problem with two possible out- 
comes. By placing Boolean variables on each node of the 
system, one can study cither equilibrium or noncquilib- 
rium properties of the thermodynamics of the Boolean 
variables and thus see how stable a given state of the 
system is. Generalizations to more complex decision 
problems can be readily accomplished by replacing the 
Boolean variables with, e.g., g-statc Potts variables [Ij. 

The entities interacting on a real network may have 
complex interactions, but models typically focus on con- 
nectivity in randomly-occupied networks or on networks 
with uniform interactions between the entities. As illus- 
trated above and as suggested in different studies 0] , 
many real networks possess both friend and foe interac- 
tions among the degrees of freedom. Because the network 
intrinsically has loops, this leads to frustration between 
the Boolean variables, quickly complicating the study of 
such systems. 

Here, we study the critical behavior of the the random- 
bond Ising model on scale-free networks. The model 
maps directly onto a friend/foe network (random bonds) 
with Boolean variables (Ising spins). Crucially, both 
ferromagnetic and antiferromagnetic interactions are al- 
lowed. Although many of the networks of experimental 
importance are dynamic and out-of-equilibrium, a thor- 
ough understanding of the equilibrium model provides 
a first step into the understanding of generic problems 
associated with networks with random interactions. 

Monte Carlo simulations of Ising spins on scale-free 
networks and complex random graphs with uniform an- 
tiferromagnetic interactions j8l4lC| have shown that a 
stable spin-glass phase exists. Similarly, studies of a 
random-field ferromagnetic Ising model [ll|, [T2| on scale- 
free graphs [l3| show that for A < 3 the spins are always 
ordered (i.e., consensus is stable to local perturbations), 
whereas for A > 3 a phase transition between a param- 
agnetic and ferromagnetic phase exists as a function of 
the random-field strength. 

Surprisingly, for the case of a pure Ising spin glass de- 
fined on a scale-free graph no detailed numerical results 
exist with most results relying on analytical approxima- 
tions and mean- field calculations |14| - [l7| . The detailed 
mean-field study by Kim et al. [l^ showed that for A < 3 
the critical temperature of the system diverges, i.e., the 
spins arc stable to arbitrary local perturbations, whereas 
for A > 3 a finitc-tcmperaturc transition from a para- 



magnetic to a spin-glass state exists. In this work we 
improve on the results by Kim et al. by expanding the 
approach of Leone et al. for ferromagnetic systems to 
spin glasses. We present analytical results backed up by 
numerical results (using large-scale Monte Carlo simula- 
tions) for both Gaussian-distributed and bimodal edge 
weights between the Ising spins. We show that when 
the probability p of positive (ferromagnetic) interactions 
and the probability 1 —p for negative (antiferromagnetic) 
interactions is the same, the system undergoes a finite- 
temperature spin-glass transition if A > 3, in agreement 
with previous results [isj . Furthermore, we show that 
spin glasses on scale-free networks seem to obey univer- 
sality. This means that, in the pure spin-glass case, the 
type of the interactions does not seem to affect the na- 
ture of the order when opinion forms. However, when 
A < 3, a spin-glass phase is stable for all temperatures. 
Finally, when both spin-glass and ferromagnetic order are 
present {p > 0.5), only spin-glass order is present. Relat- 
ing back to the social Gedankenexperiment, this would 
suggest that for certain networks the opinion of the indi- 
vidual is robust towards local perturbations and cannot 
be affected by global consensus. 

In Sec. [H] wc introduce the Hamiltonian studied, fol- 
lowed by how the networks are constructed in Sec. lIIII In 
Sec. IIVI we present analytical results and construct a A- 
p (network strength vs fraction of ferromagnetic bonds) 
phase diagram. Details about the simulations are shown 
in Sec. |Vl followed by numerical results in Sec. IVII and 
concluding remarks. 



II. MODEL 

The Hamiltonian of the Edwards- Anderson Ising spin 
glass [l^, [131 defined on a scale-free graph is given by 



N 



Kj 



(2) 



where the Ising spins Si £ {±1} lie on a scale-free graph 
with N sites and interactions 



1) 



•(3) 



In Eq. (j3|) Eij = 1 if a bond is present between spin Si and 
Sj and £ij = otherwise, n denotes the mean connectiv- 
ity of the underlying graph. The connectivity of site i, 
ki := Eij, is sampled from the scale-free distribution, 
Eq. ([T|). The bond values arc drawn from either a Gaus- 
sian distribution with zero mean and standard deviation 
unity, i.e.. 



pj(J„-) ~ exp(-j2/2) 



or a bimodal distribution defined via 



(4) 



P.Mj) = (1 - P)5{J-.j + 1) + P5{J^3 - 1) ■ (5) 

In Sec, mil wc describe in detail how the scale-free graphs 
as shown in Fig. [T]arc generated for the simulations. 
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III. GENERATING SCALE-FREE GRAPHS 

One standard approach for the generation scale- 
free networks is preferential attachment [2l|. In this 
physically-inspired growth process, a new node is added 
to the graph at each step, and the probability of attach- 
ing to previously existing nodes depends on their edge de- 
grees. This is believed to mimic the creation of scale- free 
networks in a wide variety of processes, where newcomers 
are more likely to associate with already-popular mem- 
bers of the network. The frustration that must be present 
for the disordered problem to be nontrivial requires that 
there be loops present in the network. The new nodes 
must therefore attach to multiple pre-existing nodes in 
this particular growth process. The simplest implemen- 
tation of preferential attachment, where the probability 
of attaching to a node of edge degree k is proportional to 
k produces a power law distribution of edge degrees with 
exponent A = 3. It is possible to modify the exponent, 
at least in the N ^ oo limit, by changing the function 
giving the probability of attachment [22, [23| . 

Another technique for generating scale-free networks 
is to create "annealed graphs" where an edge degree dis- 
tribution is chosen randomly according to a power law 
of exponent A and a graph is chosen randomly from the 
ensemble of all possible graphs consistent with the cho- 
sen edge degree distribution 2J]. Contrary to the name, 
once a graph is chosen, it is fixed in time for a given 
sample. This is done by assigning k stubs for each node, 
where k is drawn from the distribution pk [Eq. ([U], and 
randomly pairing the stubs. These pairings make up the 
edges in the graph. If the resulting graph is valid (in our 
case, we do not allow double edges, and only connected 
graphs arc considered), then it is accepted and may be 
used for simulation. 

In practice, wc use a slightly different approach which 
is much faster but is known to cause the selection of 
graphs to be slightly nonuniform [2^, [2^ . If, during the 
pairing process, a connection is to be made between two 
stubs corresponding to the same node, this is not allowed. 
In the method described above, all edges are removed and 
the pairing starts from the beginning. Here, we simply 
reject the pairing and move on. This is not expected to 
afl^ect our results significantly: The degree distribution 
is fixed independently of this method. In practice, the 
preferential attachment graphs are quite different than 
the annealed graphs used here, yet our tests give qual- 
itatively similar results for the two cases. The results 
presented in this paper are from simulations using an- 
nealed graphs as defined above. 

This graph-generation technique works for general de- 
gree distributions, although the acceptance rate may be 
prohibitively low for some graphs. For application to 
scale-free graphs, an upper bound is imposed on the al- 
lowed edge deg rees, /cmax — V^lV. Although we can gener- 
ate graphs with k exceeding a/ZV, the ensemble is poorly 
defined in this case: Even randomly-chosen graphs can- 
not be uncorrelated [2J, Sal- We also set a lower 



bound on the edge degree fcniin = 3. This eliminates spins 
which could be easily integrated out of the system and 
do not contribute to the frustration properties: dangling 
spins with only one attachment, and (possibly long) loops 
of spins which are not connected with any other spins. 



IV. ANALYTICAL RESULTS 

Analytical results for spin glasses on scale-free net- 
works were obtained previously jisj , and here we adapt a 
calculation for the Ising ferromagnet on a scale- free graph 
to the spin-glass case |18j . 



A. Replica approach 

We use the replica approach and at first consider the 
disorder average of the replicated partition function 
for integer powers of n 



(6) 



N 



i<j j=l 



The double angular brackets ((■■■)) denote an average 
only over the quenched interaction variables Jij and eij. 
At the moment the connectivities ki are fixed. The con- 
straints on these quantities, which are necessary to im- 
pose a scale- free degree distribution, will be introduced 
later [see Eq. (dU])] as done in Ref. [13. In Eq. (H]) A/" is a 
normalization constant and /3 = 1 /T the inverse temper- 
ature. After an appropriate continuation to small values 
of n the expression in Eq. ([5]) is related to the free energy 
F of the model via 



/3F 



lim 



iZ"')) - 1 



(7) 



Using a representation of the (5-function, the integrals in 
Eq. ([U factorize resulting in 



(^"))=E 



-K.N/2 
'JT 



N 



4=1 



X exp 



where Si is an Ising spin with n components. The order 
parameters 



(9) 



and their conjugated fields /3((t) allow one to perform 
the trace over the spin variables Si in Eq. ([5]). After 
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replacing the connectivity-dependent site averages with 
the appropriate averages over the degree distribution pj, 
we obtain 

1 ^ 

-Y^\ogY^{p[a)t ^ ^p,.logX:(/5(?))' 

i—1 a k ff 

= Pfxip)- (10) 
The partition function {{Z'"-J) then acquires the form 



\{ dp{<j)dp{cT) CXp ( ~ iV/3/tHal(p, P)) (11) 



with the trial free energy 

/3/triai(p,/5) = «^p(a)p(a) + |-/3/A(p) (12) 



where the average over the distribution pj is represented 
as (•••),/. The integral in Eq. ([TT|) can be evaluated 
with the method of steepest descents leading to a self- 
consistent equation for the order parameters p: 



k 



k-l 



(13) 



Finding stable solutions to this equation in the ti — )■ 
limit would lead to the free energy of the model at all 
temperatures. This problem remains to be solved for 
general spin-glass Hamiltonian. However, here we are 
interested only in the transition temperature of the sys- 
tem. In this case the simplest replica symmetric analysis 
is sufficient. 



B. Replica-symmetric solution 

Due to the Boolean nature of the Ising spins, the 
replica-symmetric solution Prs('?) of Eq. (jl3p only de- 
pends on the sum s = J^a '^a of the components of u. 
The parametrization 



p,^{a) ^ / dhV{h)cxp{l3hs), 



(14) 



allows one to perform the limit n — > straightforwardly 
which leads to the self-consistent equation for the local- 
field distribution 



, k-l 



k=k„ 



k f 

V{h) = V pk- / \\dhMV{h{)pj[J,) 

'— i=l 

x5(h~Y.u{h,,JM (15) 



with u defined as 

u{h, J) = iatanh[tanh(/?/i) tanh(/3J)]. (16) 

P 

As pointed out in Ref. [l^ and calculated in Ref. [2^ this 
equation can be derived within the cavity framework [soj . 
When connecting a new site to the system, one has to 
take into account the heterogeneity of the graph, which 
is reflected in the distribution [k/ K)pk on the right-hand 
side of Eq. (fT5|) . 

The (5-function is always a solution of the self- 
consistent equation. Due to the vanishing of all local 
fields, this solution intuitively corresponds to the param- 
agnetic phase. Based on physical grounds, we expect 
this solution to be unstable below a critical tempera- 
ture, which thus signals a transition to a frozen low- 
temperature phase. 

We first test the stability of this paramagnetic solution 
towards a spin-glass transition. To this end we plug in a 
distribution Vq with a small second moment (/i^)o = e on 
the right-hand side of Eq. . If the second moment of 
the resulting distribution on the left-hand side exceeds e 
the paramagnetic solution becomes unstable and the sys- 
tem undergoes a spin-glass transition. The corresponding 
critical temperature T^^ is given by the equation 



1 



E 



k{k-i) 



pk ^tanh^ 



J 



(17) 



To detect a transition towards a ferromagnetic phase a 
similar procedure using a distribution Vq with a small 
first moment (/i)o = e can be applied. This leads to the 
stability criterion to determine Tj^ 



1 



E 



k{k-i) 



pk ^tanh 



J 



J^FM 



(18) 



When lowering the temperature the system enters the 
frozen phase with the higher Tc- Within the low- 
temperature phase no conjectures can be made from this 
"paramagnetic" analysis. For a bimodal bond distribu- 
tion [Eq. ([SJ] with minimum connectivity fcmin = 3 the 
solutions of Eqs. (|17p and arc visualized in Fig. [21 
Note that both transition temperatures diverge for A < 3, 
because the second moment of the degree distribution is 
infinite. 



C. Comparison to the static model 

We briefly compare our result to previous calculations. 
In the pioneering work of Kim et at. [lB| the "static" 
model was used. In this approach a probability tt^ oc 
i^^^ p G [0,1) is assigned to each vertex i = 1...N 
to obtain a scale- free graph with degree A = 1 -I- 
and mean connectivity k. The results derived for the 
phase-boundaries in Ref. [l^ rely on a truncation scheme 
which conflnes the number of the order parameters to the 
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FIG. 2: (Color online) Analytical phase diagram for fraction 
of ferromagnetic bonds p and exponent A. For A < 3, both 
T^^'^ and T^*^ diverge, but there is a crossover in the rates 
at which they diverge. In this case, in the blue-white region 
(top) the system is ferromagnetic at all finite temperatures, 
while in the red-yellow region (bottom) , the system is a spin 
glass at all finite temperatures. Note that for A — >■ 2 the sys- 
tem becomes increasingly stable against local (temperature) 
perturbations. 



"most important" ones [31|. At the replica-symmetric 
level, this approximation is equivalent to the assumption 
of Gaussian local fields, which can be avoided, using a 
variant of the approach described above. 

Starting from Eq. (10) in Rcf. [IE! we proceed using the 
order parameters 



(19) 



and their conjugated variables p{a). The functions which 
still depend on the vertex-weights tt^ can be transformed 
to integrals in the thermodynamic limit according to 



p{k)g{k). 



In the last step we performed a substitution which di- 
rectly leads to the scale-free distribution relevant for the 
static model 



P(fc) 



, A-l 



A-1 



A 



A-l 



K. (20) 



The partition function ((.^")) of the static model reduces 
to a saddle-point integral [see Eq. pT|) ] with a scale-free 
dependent part in the the trial free energy 



poo 



dk p{k) log 



Using the saddle-point equations and the replica- 
symmetric Ansatz, we obtain a self-consistent equation 
for the distribution of local fields 



r{h) = rdkp{k)^Y.^'^ 



K ^ — ' m! 

m=0 



Y[dhM (21) 



1=1 



This equation is a generalization of the replica-symmetric 
equations for the Viana-Bray model [Sll, 132] where the 
mean connectivity k is sampled from the distribution 
{k/ K)p{k). Moreover, the equation is similar to the self- 
consistent equation [Eq. ((T5|) ]. where the connectivities 
are sampled from the distribution {k/ K)pk. This last ap- 
proach is a generalization of the fixed-connectivity model 
to arbitrary degree distributions, and we prefer it due to 
its generality. The paramagnetic solution V{-) = 5{-) is 
a solution of Eq. (PT|) which becomes unstable towards a 
spin-glass transition at the critical temperature 



1 



K ( tanh — ^ 



dk p{k)- 

°(A-2)^ 
(A-l)(A-3)• 



(22) 



The transition to a ferromagnetic phase occurs at the 
critical temperature 



1 = K ( tanh , 

\ V T™ 



(A -2)2 



(A-l)(A-3) 



(23) 



In the last line of Eq. we performed the integrals 
which are convergent for A > 3 only. This representation 
coincides with the predictions for the critical tempera- 
ture in Ref. ITsl . Here the results were obtained without 
resorting to a truncation scheme. Predictions for the de 
Almeida-Thouless line [33| relying on Gaussian approxi- 
mations |34| can lead to results which were shown to be 
wrong [33, l3a|. The investigation of the phase bound- 
aries in the presence of an external magnetic field seems 
promising in this approach, and will be reported else- 
where. 



V. NUMERICAL DETAILS 

We validate the aforementioned analytical results using 
Monte Carlo simulations. In particular, the numerical 
results show the strength of the different corrections to 
scaling depending on the choice of the exponent A of the 
scale-free network. 



A. Observables 

To determine the location of both the ferromagnetic 
and spin-glass phase transitions we measure first the 
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Binder cumulant [33 defined via 



2 I (((O^))) 



(24) 



In Eq. (P^ . (((• • •))) represents an average over the disor- 
der via (• • • ) J, the random graphs via {■ ■ ■)k, and Monte 
Carlo time, i.e., (•••)t- Furthermore, O is either the 
magnetization m in the ferromagnetic sector defined via 



(25) 



or the spin-glass overlap q in the spin-glass sector given 
by 



1 ^ 



(26) 



In Eq. (|26l). "a" and "/?" represent two copies of the 
system with the same disorder. The Binder ratio is di- 
mensionless and thus has the simple scaling form 



g = G(L^/''[T-Tc 



(27) 



where Tc represents the transition temperature. Note 
that the two-point finite-size correlation length [38l - l4l| 
typically used to pinpoint transitions in glassy systems 
cannot be used here because the spins are placed on a 
lattice that has no geometry. 

The calculation of the Binder ratio g allows one to 
determine Tc and the critical exponent v for both the 
spin-glass and ferromagnetic sector. However, to fully 
characterize the critical behavior of the model a second 
critical exponent has to be computed [H. We have also 
computed the susceptibility x given by 



where the scaling behavior 



(28) 



(29) 



yields the exponent 77. However, fluctuations in the data 
are huge and thus the determination of the critical expo- 
nent 77 was not possible. 



B. Equilibration and simulation parameters 

The simulations are done using the parallel tempering 
Monte Carlo method [13]. For p = 0.5 we first simulated 
the system with Gaussian disorder to obtain an idea of 
the equilibration behavior of the system. Furthermore, 
in the Gaussian case we can perform a rigorous equili- 
bration test [H, ll^l where the energy is compared to an 
expression derived from the link overlap, sec Eq. pop . 
The data for both quantities approach a limiting value 
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FIG. 3: (Color online) Example equilibration according to 
Eq. ((301) for ^ = 8192 spins at T = 1.027 (lowest temperature 
simulated) for A — 4.5 and Gaussian bonds between the spins. 
Once the data for the energy U and the energy computed from 
q4 [f/(g4)] agree, the system is in thermal equilibrium. This 
can be seen also with data for that also are independent of 
Monte Carlo time (Note that the data are shifted by 1.2 for 
better viewing). 



from opposite directions. Once they agree, and other 
observables are independent of Monte Carlo time, the 
system is in thermal equilibrium. 

When the values of the interactions between the spins 
arc drawn from a Gaussian distribution, one can use the 
following expression to test if the Monte Carlo simula- 
tions have converged to thermal equilibrium: 



Uiq4) 



1 

T 



N 



(1 - 94) 



where 



94 



„a a B B 

i^i Sj Sj Sj 



(30) 



(31) 



and Nb is the number of nonzero bonds of a given sam- 
ple. Note that the expression in Eq. pop is valid for the 
spin-glass sector, however, it represents a conservative 
bound for the ferromagnetic sector [4^. Furthermore, 
Nf, is inside the disorder average because the number of 
bonds fluctuates from sample to sample. Sample data are 
shown in Fig. [S] Once both sides of Eq. ((30)) agree, the 
data for the squared order parameter (shifted for better 
viewing) is also in thermal equilibrium. 

For the bimodal disorder distribution and the ferro- 
magnetic sector the aforementioned equilibration test 
cannot be used. In this case we preform a logarithmic 
binning of all observables. Once the data for the last 
three bins agree within error bars we deem the system 
to be in thermal equilibrium. Simulation parameters are 
shown in Table [H 



TABLE I: Parameters of the simulation: For each exponent A 
and fraction of ferromagnetic bonds p (Note: Gaussian disor- 
der is marked with 'Gauss') we compute A'sa disorder/network 
instances. Ns^ = 2*" is the total number of Monte Carlo 
sweeps for each of the 2Nt replicas for a single instance, Tmin 
[Tinax] is the lowest [highest] temperature simulated, and Nt 
is the number of temperatures used in the parallel tempering 
method for each system size A*'. 
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VI. NUMERICAL RESULTS 



We first study Gaussian disorder where we have a 
strong equilibration test to ensure that the data are in 
thermal equilibrium. Furthermore, by scaling the data 
of the Binder ratio we can study the critical exponent v 
and subsequently compare to results for bimodal disorder 
to test for universality. Corrections to scaling are very 
large for this model despite the large number of samples 
studied. Therefore, a finite-size scaling analysis of the 
susceptibility has not been done, which is why any claims 
for universal behavior are indicative only (two exponents 
are needed for a strong statement). 
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FIG. 4: (Color online) Binder ratio Qq for the SG sector and 
A = 3.0 with Gaussian disorder for different system sizes N . 
The crossings between increasing system-size pairs diverge 
with increasing system size suggesting that T^''^ — ;> 00, in 
agreement with the analytical predictions. 



A. Gaussian disorder 

Figure 2] shows data for the Binder ratio gq for Gaus- 
sian disorder and A = 3.0, right at the onset (see Fig. [5]) 
where the critical temperature for the spin-glass (SG) 
sector starts to diverge. The crossing temperatures be- 
tween lines for N/2N pairs grow with the system size 
in agreement with the analytic calculations. To prevent 
rpSG -j-j^g gQ sector to diverge, the bonds would have 
to be re-scaled. Furthermore, there is no transition in 
the ferromagnetic sector (not shown), in agreement with 
the analytical calculations. 

In Fig. Owe show data for A = 3.5. In agreement with 
the analytical predictions Tf^ is finite, albeit with huge 
corrections to scaling. We estimate T^'^ = 2.31(3). Note 
that this estimate is computed via a finite-size scaling of 
the data and only takes statistical fluctuations into ac- 
count. We have no control over finite-size corrections. A 
crude extrapolation suggests that the critical tempera- 
ture will likely be larger than the quoted value which we 
treat as a lower bound. This is a generic problem for the 
estimates of the critical temperature made in the spin- 
glass sector that decreases drastically for increasing A. 
Again, there is no transition in the ferromagnetic sector. 

Finally, Fig. [5] shows data for the Binder ratio for A = 
4.5. The data cross at T^'^ ^ 1.39(1). For all values of 
A studied where Tc < 00 the estimates for the critical 
exponent i' are listed in Table [TTl Although statistical 
error bars are large (up to ^ 16% of the estimate in the 
Gaussian case), the results suggest that for the values 
of A studied the systems might share the same critical 
exponent v. This, in turn, suggests that the different 
values of A might share the same universality class for 
A > 3.0. 
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FIG. 5: (Color online) Binder ratio Qq for the SG sector and 
A = 3.5 with Gaussian disorder for different system sizes TV. 
The data show strong corrections to scaling. Using a finite- 
size scaling analysis we estimate T^*"" = 2.31(3). In this and 
all subsequent figures, the width of the shaded region around 
the critical temperature corresponds to the statistical uncer- 
tainty for the estimate of the critical temperature using a 
combination of a finite-size scaling analysis and a bootstrap 
method. 



10" 



10- 



10- 



H»H JV 


= 512 


■ JV 


= 1024 


^ JV 


= 2048 


M»M JV 


= 4096 


^•^ JV 


= 8192 




A = 4.5 

= 1.39(1) 
Gaussian — SG sector 



1.1 1.2 1.3 1.4 

T 



1.5 1.6 



FIG. 6: (Color online) Binder ratio Qq for the SG sector 
and A = 4.5 with Gaussian disorder for different system 
sizes A'^. Corrections to scaling are small and we estimate 
= 1.39(1). 



B. Bimodal disorder 

Figure [7] shows data for the Binder ratio gq for bimodal 
disorder and A = 3.0 and p = 0.50, right at the onset (see 
Fig. [5]) where the critical temperature for the spin-glass 
(SG) sector starts to diverge. As for the Gaussian case 
presented in Sec. IVI Al crossing points between lines for 
N /2N pairs grow with the system size in agreement with 
the analytic calculations, i.e., T^'^ — )• oo. Furthermore, 
for the ferromagnetic sector there is no transition (not 
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FIG. 7: (Color online) Binder ratio Qq for the SG sector and 
A = 3.0 with bimodal disorder for different system sizes A'' 
and p = 0.50. As for the Gaussian case (Fig.|4]), the crossings 
between increasing system-size pairs diverge with increasing 
system size suggesting that T^*"" — >■ oo, in agreement with the 
analytical predictions. 
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FIG. 8: (Color online) Binder ratio Qq for the SG sector and 
A = 3.5 with bimodal disorder for different system sizes A' 
and p = 0.50. As for the Gaussian case (Fig. [S]), the data 
show strong corrections to scaling. Using a finite-size scaling 



analysis we estimate 



2.55(8). 



shown) . 

In Fig. [5] we show data for A = 3.5 and p = 0.50. 
In agreement with the analytical predictions T^'^ is fi- 
nite, albeit with large corrections to scaling. We esti- 
mate T^'^ = 2.55(8). Again, there is no ferromagnetic 
transition (not shown). 

We now study in detail A = 4.5 for different concen- 
trations of ferromagnetic bonds p. Figure ^ shows data 
for the Binder ratio for A = 4.5 and p = 0.50 for both 
spin-glass (SG) and ferromagnetic (FM) sectors. For the 
spin-glass sector (top panel), the data cross cleanly at 
T^SG ^ 1.70(1). However, for the FM sector (bottom 
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FIG. 10: (Color online) Top panel: Binder ratio gq for the 
SG sector and A = 4.5 with bimodal disorder (p — 0.70) for 
different system sizes N. We estimate rj^'^ = 1.69(1). Bottom 
panel: Binder ratio Qm for the FM sector. The shaded area 
in the bottom panel corresponds to T^^ — 1.70 ± 0.01 (top in the bottom panel corresponds to T^'~' — 1.69 ± 0.01 (top 



FIG. 9: (Color online) Top panel: Binder ratio gq for the 
SG sector and A = 4.5 with bimodal disorder (p — 0.50) for 
different system sizes A'^. We estimate T^'^ = 1.70(1). Bottom 
panel: Binder ratio Qm for the FM sector. The shaded area 

^SG 



panel). The data show no sign of ferromagnetic order: The 
data do not cross and decrease for increasing system size A'". 



panel). The data show no sign of ferromagnetic order: The 
data do not cross and decrease for increasing system size A^. 



panel) the data strongly suggest that Tj'^ is not de- 
fined, in agreement with the analytic predictions. For 
A = 4.5 the phase boundary between the FM and the 
SG sector lies somewhere between p = 0.7 and 0.85, see 
Fig. [21 Therefore, we expect that for p = 0.7 we only 
have SG order, whereas for p = 0.85 the system orders 
ferromagnctically. 

The case for p — 0.70 is shown in Fig. [TUl Like 
for p — 0.50, the data for the SG sector show a clear 
transition, whereas the data for the FM sector suggest 
that there is no transition. Note that the estimate for 
= 1.69(1) agrees within error bars with the estimate 
for p = 0.50, suggesting that the phase diagram does not 
depend on p in the SG regime, in agreement with the 
analytical calculations. For p = 0.85 there is no SG or- 
der (not shown, in agreement with the analytical results). 
Moreover, T™ = 2.428(8), see Fig.HIl Interestingly, for 
all values of p where < T^^*^ < oo the critical exponent 



V agrees within error bars, suggesting universal behavior. 
However, the error bars are very large and a different sce- 
nario might emerge from better data. Note also, that for 
the individual values of A studied, the critical exponent v 
agrees for both bimodal and Gaussian disorder, see Table 

El 



C. Universality? 

To determine if two systems share the same universal- 
ity class two independent critical exponents need to be 
computed. However, fluctuations in the data are large 
and therefore estimating the critical exponent 77 from 
a finite-size scaling analysis of the susceptibility is dif- 
ficult. The estimates for the critical exponent v agree 
within error bars for all values of A studied, independent 
of the used disorder (Gaussian versus bimodal). This 
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FIG. 11: (Color online) Binder ratio gm for the FM sector. 
The shaded area corresponds to T™ = 2.428(8). In this case 
we are deep in the FM sector, the SG sector shows no sign of 
a transition (not shown). 



TABLE II: Critical parameters Tc and f for the spin glass 
(SG) and ferromagnetic (FM) sectors computed using a finite- 
size scaling analysis of the data for A'^ > 2048. 
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would, at first sight, give a strong indication of uni- 
versal behavior. To verify this result, we compare a 
second universal quantit y, t he value of the dimension- 
less Binder ratio at [4^. For A = 3.5, 5(Tc,A 
3.5, Gauss) — 0.018(3), which agrees within error bars 
with .g(Tc,A = 3.5, p = 0.5) = 0.012(25). For A = 4.5, 
g{Tc,X = 4.5, Gauss) = 0.132(12), which agrees within 
error bars with g{Tc,\ = 4.5, p = 0.5) = 0.129(12). In 
the bimodal case, we find also that g{Tc,X = 4.5, p = 



0.5) = 0.129(12) agrees within error bars with g(Tc,X = 
4.5, p = 0.7) — 0.136(12), thus suggesting the the same 
universality class might be shared below the spin-glass- 
to-ferromagnet phase boundary for all values of p. This 
suggests that for a given value of A the networks with 
both Gaussian and bimodal disorder might share the 
same universality class. 



VII. CONCLUSIONS 

We have studied Boolean (Ising) variables on a scale- 
free graph with competing interactions. Our analytical 
and numerical results show that for A < 3 the critical 
temperature diverges with the system size. For larger val- 
ues of A the system undergoes a finite-temperature tran- 
sition between a spin-glass and a paramagnetic phase. 
The robustness of both the ferromagnetic and spin-glass 
phases suggest that Boolean decision problems on scale- 
free networks are quite stable to local (temperature) per- 
turbations. For the case with bimodal disorder, we show 
that for a large enough fraction of ferromagnetic bonds 
the system orders ferromagnetically at finite tempera- 
tures. 

Real networks typically have exponents A < 3. It 
would be interesting to study such a network in the fu- 
ture. Furthermore, the effect of "global biases" (field 
terms) will also be studied. Finally, opinion formation is 
an intrinsically nonequilibrium process, a problem that 
will also be studied in the future. 
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